Measurement of Day and Night Neutrino Energy Spectra at SNO and 
Constraints on Neutrino Mixing Parameters 
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The Sudbury Neutrino Observatory (SNO) has measured day and night solar neutrino energy 
spectra and rates. For charged current events, assuming an undistorted *B spectrum, the night 
minus day rate is 14.0% ± &.'i%^W% of the average rate. If the total flux of active neutrinos is 
additionally constrained to have no asymmetry, the Ve asymmetry is found to be 7.0% ± 4.9%1}'2%. 
A global solar neutrino analysis in terms of matter-enhanced oscillations of two active flavors strongly 
favors the Large Mixing Angle (LMA) solution. 
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PACS numbers: 26.65.+t, 14.60.Pq, 13.15.+g, 95.85.Ry 



The Sudbury Neutrino Observatory (SNO) has pro- 
vided strong evidence that neutrinos change flavor as 
they propagate from the core of the Sun ||, inde- 
pendently of solar model flux predictions. This flavor 
conversion can be explained by neutrino oscillation mod- 
els based on flavor mixing. For some values of the mixing 
parameters, spectral distortions and a measurable depen- 
dence on solar zenith angle are expected |§, The 
latter might be caused by interaction with the matter of 
the Earth (the MSW effect) and would depend not only 
on oscillation parameters and neutrino energy, but also 
on the path length and e~ density through the Earth. 
This Letter presents SNO's first measurements of day 
and night neutrino energy spectra, and constraints on 
allowed neutrino mixing parameters. 

The data reported here were recorded between Novem- 
ber 2, 1999 and May 28, 2001 UTC. The total livetimes 
for day and night are 128.5 and 177.9 days, respectively. 
The time-averaged inverse-square distance to the Sun 
was 1.0002 (day) and 1.0117 (night). During 
the development of this analysis, the data were parti- 
tioned into two sets of approximately equal livetime (split 
at July 1, 2000), each having substantial day and night 
components. Analysis procedures were refined during the 
analysis of Set 1 and fixed before Set 2 was analysed. The 
latter thus served as an unbiased test. Unless otherwise 
stated, the analysis presented in this paper is for the 
combined data set. 

The data reduction in Q was used here. For each 
event, the number, pattern, and timing of the hit photo- 
multiplier tubes (PMTs) were used to reconstruct effec- 
tive recoil electron kinetic energy T^g, radial position R, 
and scattering angle 6q with respect to the Sun-Earth 
direction. The charged current (CC), elastic scattering 
(ES) and neutral current (NC) reactions each have char- 
acteristic probability density functions (pdfs) of Teg, R, 
and 9q . A maximum likelihood fit of the pdfs to the data 
determined the fiux from each of these reactions. 

The measured night and day fluxes (fiN and (fin were 
used to form the asymmetry ratio for each reaction: A = 
2(0Ar — (f)D)/{4'N + (po)- The CC interaction is sensitive 
only to i^e- The NC interaction is equally sensitive to all 
active neutrino flavors, so active-only neutrino models 
predict Anc = 1^. The same models allow Acc 7^ 0. 
The ES reaction has additional contributions from v^r 
leading to a reduction in its sensitivity to asymmetries. 

SNO used calibration sources to constrain varia- 
tions in detector response pi that can lead to day-night 
asymmetries. A ^^N source g, which produces 6.1-MeV 
gamma rays, revealed a 1.3% per year drift in the energy 
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FIG. 1: Various event classes used to determine systematic 
differences between day and night measurements. Also shown 
are measured asymmetries on the CC flux, and on the electron 
neutrino flux derived from the CC, ES, and NC rates when 
the total neutrino flux is constrained to have zero asymmetry. 



scale. Due to seasonal variation in day and night livetime, 
this drift can create an artificial asymmetry. The analy- 
sis corrected for this drift and a systematic uncertainty 
was assigned using worst-case drift models. Gamma rays 
from the ^^N source were also used to constrain direc- 
tional dependences in SNO's response. 

A set of signals that are continuously present in the 
detector was used to probe possible diurnal variations in 
detector response. The detector was triggered at 5 Hz 
with a pulser, verifying livetime accounting. Muons pro- 
vide an almost constant signal and, through interactions 
with D2O, produce secondary neutrons. After applying 
a cut to remove bursts with high neutron multiplicity, 
these muon-induced neutrons were used to limit temporal 
variations in detector response. A more sensitive study 
focused on a solitary point of high background radioac- 
tivity, or "hot spot" , on the upper hemisphere of the SNO 
acrylic vessel, apparently introduced during construction. 
Its event rate was stable and sufficient to make an excel- 
lent test of diurnal variations. It also provides a sensitive 
test for changes in reconstruction. A limit of 3.5% on 
the hot spot rate asymmetry was determined, which be- 
cause of its steeply falling energy spectrum constrained 
the day and night energy scales to be the same within 
0.3%. An east/west division of the neutrino data based 
on the Sun's position should show no rate variations from 
matter effects. As expected, the CC rates for east and 
west data were consistent. The rate asymmetries for each 
test are shown in Fig. |l|. 

Backgrounds were subtracted separately for day and 
night as part of the signal extraction. The results were 
normalized for an Earth-Sun distance of 1 AU, yielding 
the results in Table |. Day and night fiuxes are given 
separately for data Sets 1 and 2, and for the combined 
data. A consistency test of the six measured fiuxes 
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TABLE L The results of signal extraction, assuming an undistorted *B spectrum. The systematic uncertainties (combined 
set) include a component that cancels in the formation of the A. Except for the dimensionless A, the units are 10® cm~^ s~^. 
Flux values have been rounded, but the asymmetries were calculated with full precision. 



signal 


Set 1 

015 (t>N 


Set 2 

015 0JV 


Combined 

0D 0JV 


A{%) 


CC 
ES 
NC 


1.53 ±0.12 1.95 ±0.10 
2.91 ±0.52 1.59 ±0.38 
7.09 ±0.97 3.95 ±0.75 


1.69 ±0.12 1.77 ±0.11 
2.35 ±0.51 2.88 ±0.47 
4.56 ±0.89 5.33 ±0.84 


1.62 ±0.08 ±0.08 1.87 ±0.07 ±0.10 
2.64 ±0.37 ±0.12 2.22 ± 0.30 ± 0.12 
5.69 ± 0.66 ± 0.44 4.63 ± 0.57 ± 0.44 


±14.0 ± G.Stil 
-17.4± 19.5t2;2 
-20.4± 16.9l2l 




Kinetic energy (MeV) Aq (%) 



FIG. 2: (a) Energy spectra for day and night. All signals 
and backgrounds contribute. The final bin extends from 13.0 
to 20.0 MeV. (b) Difference, night - day, between the spectra. 
The day rate was 9.23 ± 0.27 events/day, and the night rate 
was 9.79 ± 0.24 events/day. 



FIG. 3: Joint probability contours for Atot and Ae- The 
points indicate the results when Atot is allowed to float and 
when it is constrained to zero. The diagonal band indicates 
the 68% joint contour for the Super-K Aes measurement. 



changing variables to (jjcc = 4>e, 4>nc = 4>tot 



between Sets 1 and 2 yielded a chance probability of 
8%. A similar test done directly on the three asym- 
metry parameters gave a chance probability of 2%. No 
systematic has been identified, in either signal or back- 
ground regions, that would suggest that the differences 
between Set 1 and Set 2 are other than a statistical fluc- 
tuation. For the combined analysis, Acc is ±2.2(t from 
zero, while Aes and Anc are ~0.9a and —1.20- from 
zero, respectively. Note that Acc and Anc are strongly 
statistically anti-correlated {p = —0.518), while Acc and 
■^ES {p — —0.161) and Aes and Anc {p = —0.106) are 
moderatedly anti-correlated. Table |l| gives the system- 
atic uncertainties on the asymmetry parameters. The 
day and night energy spectra for all accepted events are 
shown in Fig. 0. 



Table III (a) shows the results for Ae derived from the 
CC day and night rate measurements, i.e., Ae = Acc- 
The day and night flavor contents were then extracted by 



and 



OES 



where e = 1/6.48 is the ratio 



of the average ES cross sections above 5 MeV for 
and i^e- Table [II (b) shows the asymmetries of (jje and 
(j)tot with this additional constraint from the ES rate 
measurements. This analysis allowed for an asymmetry 
in the total flux of neutrinos (non-zero Atot), with 
the measurements of Ae and Atot having a strong anti- 
correlation. Fig. II shows the Ae vs. Atot joint prob- 
ability contours. Forcing Atot = 0, as predicted by 
active-only models, yielded the result in Table III (c) of 



Ae = 7.0% ±4.9% (stat.)li;^% (sys.) 

The Super-Kamiokande (SK) collaboration measured 
Aes{SK) = 3.3% ±2.2% istat.)t\i% (sys.) The 
ES measurement includes a neutral current component, 
which reduces the asymmetry for this reaction relative 
to Ae ||ll|. Aes{SK) may be converted to an equivalent 
electron flavor asymmetry using the total neutrino flux 
measured by SNO, yielding Ae{SK) (Table R (d)). This 
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TABLE II: Effect of systematic uncertainties on A (%). 
For presentation, uncertainties have been symmetrized and 
rounded. 



Systematic 


SAcc 


SAes 


5.4 JVC 


Long-term energy scale drift 


0.4 


0.5 


0.2 


Diurnal energy scale variation 


1.2 


0.7 


1.6 


Directional energy scale var. 


0.2 


1.4 


0.3 


Diurnal energy resolution var. 


0.1 


0.1 


0.3 


Directional energy resolution var. 


0.0 


0.1 


0.0 


Diurnal vertex shift var. 


0.5 


0.6 


0.7 


Directional vertex shift var. 


0.0 


1.1 


0.1 


Diurnal vertex resolution var. 


0.2 


0.7 


0.5 


Directional angular recon. var. 


0.0 


0.1 


0.1 


PMT /3-7 background 


0.0 


0.2 


0.5 


AV+H2O /3-7 bkgd. 


0.0 


0.6 


0.2 


D2O /3-7, neutrons bkgd. 


0.1 


0.4 


1.2 


External neutrons bkgd. 


0.0 


0.2 


0.4 


Cut acceptance 


0.5 


0.5 


0.5 


Total 


1.5 


2.4 


2.4 



TABLE III: Measurement of the (jf>e and (j>tot asymmetry for 
various constraints. All analyses assume an undistorted *B 
spectrum. 



Constraints 


Asymmetry (%) 


a) no additional constraint 


Acc = 14.0±6.3li-^ 
Anc = -20.4± 16.9l^;^ 
(see text for correlations) 


b) (jiES = (1 — e)<?!'e + e0tot 


Ae = 12.8±6.2+i;^; 
Atot = -24.2± 16.ll^;^ 
correlation = -0.602 


c) (j)ES = (1 - (;)0e + e0tot 
Atot = 


Ae = 7.0 ± 4.9li;2 


d) (f)ES = (1 — <^)4>e + ^4>tot 
Atot = 

Aes{SK) = 3.3% ± 2.2%tYi% 


Ae{SK) = 5.3±3.7+^-;; 
(derived from SK .4es 
and SNO total *B flux) 



value is in good agreement with SNO's direct measure- 
ment of Ae, as seen in Fig. ||. 

SNO's day and night energy spectra (Fig. ||) have also 
been used to produce MSW exclusion plots and limits 
on neutrino flavor mixing parameters. MSW oscillation 
models [|l2| between two active flavors were fit to the 
data. For simplicity, only the energy spectra were used 
in the fit, and the radial R and direction cos^© informa- 
tion was omitted. This procedure preserves most of the 
ability to discriminate between oscillation solutions. A 
model was constructed for the expected number of counts 
in each energy bin by combining the neutrino spectrum 
p^ , the survival probability, and the cross sections iQ 
with SNO's response functions S. For this analysis, the 



dominant systematics are those for the combined fluxes, 
as detailed in Refs. IT] and and not the diurnal sys- 
tematics of Table ||. 

There are 3 free parameters in the fit: the total fiux 
(pB, the difference Am^ between the squared masses of 
the two neutrino mass eigenstates, and the mixing angle 
6. The flux of higher energy neutrinos from the solar 
hep reaction was fixed at 9.3 x 10"^ cm~^ s~^ Con- 
tours were generated in Am^ and tan^ 9 for Ax^(c.L) — 
4.61 (90%), 5.99 (95%), 9.21 (99%), and 11.83 (99.73%). 
Fig. H(a) shows allowed mixing parameter regions us- 
ing only SNO data with no additional experimental 
constraints or inputs from solar models. By includ- 
ing fiux information from the CI ||l6|| and Ga experi- 
ments [0 g g |0[ |l), the day and night spectra 
from the SK experiment [p^ , along with solar model 
predictions for the more robust pp, pep and ^Be neu- 
trino fiuxes ||l5|, the contours shown in Fig. ^(b) were 
produced. This global analysis strongly favors the Large 
Mixing Angle (LMA) region (see Table IV), and tan^ 6* 
values < 1 . While the absolute chi-squared per degree of 
freedom is not particularly large for the LOW solution, 
the difference between chi-squared values still refiects the 
extent to which one region of MSW parameter space is 
favored compared to another. Repeating the global anal- 
ysis using the total SNO energy spectrum instead of sep- 
arate day and night spectra gives nearly identical results. 

In summary, SNO has measured the day-night asym- 
metries of the CC, NC, and ES reaction rates. From 
these results the first direct measurements of the day- 
night asymmetries in the i^e flux and the total v fiux 
from the Sun have been deduced. A global fit to SNO's 
day and night energy spectra and data from other solar 
neutrino experiments strongly favors the LMA solution 
in a 2-fiavor MSW neutrino oscillation analysis. 
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US: Dept. of Energy; UK: PPARC. We thank the SNO 
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TABLE IV: Best fit points in the MSW plane for global MSW 
analysis using all solar neutrino data. (j)B is the best-fit *B 
flux for each point, and has units of 10® cm~^ s~^. Am^ has 
units of eV^. Ae is the predicted asymmetry for each point. 



Region 


xLn/dof 


<t>B 


A,{%) 


Am^ 


tan^ e 


c.L(%) 


LMA 


57.0/72 


5.86 


6.4 


5.0 X 10-" 


0.34 




LOW 


67.7/72 


4.95 


5.9 


1.3 X 10"^ 


0.55 


99.5 
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FIG. 4: Allowed regions of the MSW plane determined by a fit to (a) SNO day and night energy spectra and (b) with 
additional experimental and solar model data. The star indicates the best fit. See text for details. 
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